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ABSTRACT 



If predictions for species extinctions hold, then the 'tree of life' today may be 
quite different to that in (say) 100 years. We describe a technique to quantify 
how much each species is likely to contribute to future biodiversity, as measured 
by its expected contribution to phylogenetic diversity. Our approach considers 
all possible scenarios for the set of species that will be extant at some future time, 
and weights them according to their likelihood under an independent (but not 
identical) distribution on species extinctions. Although the number of extinction 
scenarios can typically be very large, we show that there is a simple algorithm that 
will quickly compute this index. The method is implemented and applied to the 
prosimian primates as a test case, and the associated species ranking is compared 
to a related measure (the 'Shapley index'). We describe indices for rooted and 
unrooted trees, and a modification that also includes the focal taxon's probability 
of extinction, making it directly comparable to some new conservation metrics. 

Keywords: phylogenetic diversity, extinction, biodiversity conservation, Shapley index 

Short title: Taxon-specific indices of expected future biodiversity 



- 3- 



1. Introduction 

Within a given taxonomic group, individual biological species are generally considered 
to be of equal or near-equal biodiversity value. So, for instance, areas with a greater number 
of species are more valuable than those with fewer (Myers et al. 2000). When wild species 
are ranked by value, this is usually based on their threat of extinction (see, e.g. SARA 
2002). However, as pointed out by Cousins (1991), species are discovered and identified 
because they are different from other species, which suggests that they may differ in value. 
In the context of conservation, Avise (2005) has highlighted five different currencies for 
valuing species: rarity, distribution, ecology, charisma, and phylogeny. Here, we consider 
the value of a species based on its position in a phylogeny. A phylogeny is the directional, 
acyclic graph depicting relationships between leaves (species), which we define formally 
in the next section. A phylogeny generally has a root (which assigns direction) and edge 
weights that can represent unique feature diversity (e.g. as measured by evolutionary time 
or genetic distance). Species can be defined by the features they possess, and one measure 
of their worth is their expected contribution of unique features. In this way, we can use 
a phylogeny to assign a measure of evolutionary value to a species based on its expected 
contribution of unique features. Because of the highly imbalanced shape of the Tree of Life, 
some species in a phylogeny will have far fewer close relatives than others in that phylogeny 
(Mooers and Heard 1997), and these more distantly-related species will be expected to 
contribute more unique features (Faith 1992). 

Phylogenetic measures of conservation value have a long pedigree (see, e.g. Alschul 
and Lipman 1990; May 1990) and have begun to be explored in some detail (Haake et al. 
2005; Hartmann and Steel 2007; Pavoine et al. 2005a, 20056; Redding and Mooers 2006). 
So, for example, Pavoine and colleagues presented one new measure of originality, a set of 
sampling weights such that the expected pairwise distance on the tree is maximized. Haake 
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and colleagues extended the 'Shapley value' (Shapley 1953) from co-operative game theory 
to the conservation setting to calculate the average distance of a focal species to all possible 
subsets of taxa. For both measures, more original species are those expected to contribute 
more to the resulting sets. Yet another measure that uniquely apportions the tree to its tips 
(Isaac et al. 2007) and which is the focus of a new international conservation initiative (the 
EDGE initiative, Zoological Society of London) scales almost perfectly with the Shapley 
value (unpublished results). 

One question with these measures concerns the sets that individual species are asked 
to complement. For instance, given known extinction probabilities for species, some future 
sets are much more likely than others and so some species will be more valuable because 
their close relatives are less likely to be included in future sets. Here we formalize this idea 
to extend the Shapley value of a species to include pre-assigned extinction probabilities. We 
then compare our measure with the original Shapley value using the prosimian primates as 
a test case. 

Definition Let T be a rooted or unrooted phylogenetic tree with leaf set X, together 
with an assignment of positive lengths to the edges (branches) of T. We let /(e) denote the 
length of edge e, and let E(T) denote the set of edges of T. For a subset S of X, let PD(S) 
denote the phylogenetic diversity of S defined as follows. If T is unrooted then PD(S) is 
the sum of the lengths of the edges (branches) of T in the minimal subtree that connects 
S. If T is rooted, then PD(S) is the sum of the lengths of the edges of T in the minimal 
subtree of T that connects S and the root of the tree. Figure 1 illustrates these concepts, 
and includes values at the tips that we will use in the next section. Note that although the 
branch lengths in this example are clock-like, this assumption is not required in any of the 
results we describe. 
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Fig. 1. — (a) A small rooted tree with edge lengths (of 2 units for the terminal edges incident 
with C,D,E, and 1 unit for the other five edges). Each tip j has an associated extinction 
probability P(ext) = tj. (b) For a subset S = {B, D} of taxa that are extant at some future 
time, the phylogenetic diversity score PD(S) is the sum of the lengths of the edges indicated 
in bold. The dashed edges lead to extinct taxa. 



The quantity A PD (S,i) measures how much phylogenetic diversity i contributes to the tree 
that one obtains from T once species not in S have been pruned out (for example if they 
go extinct). Alternatively, Ap^(S, i) is the marginal increase in phylogenetic diversity of S 
if i is added. 



2. 



The HED index 



For a leaf i E X, and a subset S C X — {i} let 



A PD (S, i) := PD(S U {i}) - PD(S). 



Now, suppose that each species has an associated extinction probability P(ext) 
(which may vary from species to species) — for example, this may be the probability that 
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the species is extinct in (say) 100 years from now (either globally, or in some specified 
community). We will denote this P(ext) value for species j by 6j. In this paper we consider 
the simplest model which assumes that the extinction of each species in X comprise 
independent events. Given i e X, let Si denote the random subset of species in X — {i} 
which survive (i.e. do not go extinct). 

By the independence assumption we have: 

m = S] = H(l-e ] )x J] ej. 

j£S jeX ~{i}-S 

For i e X, let ipi denote the expected value of Apo(Si,i). That is, 

il>i=E[A PD (Si,i)]= m^S]A PD (S,i). (1) 

SCX-{i} 

We call tpi the heightened evolutionary distinctiveness of species i, and the function % \— > ifti 
the heightened evolutionary distinctiveness (HED) index for T. Notice that if all the species 
in X — {i} were guaranteed to survive, then ipi would be just the length of the pendant edge 
incident with leaf i, however random extinctions mean that ipi will tend to be increased 
('heightened') over this pendant edge length. 

A related but different index, based on the Shapley value in co-operative game theory, 
has recently been described by Haake et al. (2005). This index, denoted here as ip sh can be 
defined (for unrooted trees) as follows: For i E X, 

^=\Y\ £ ( m ~ l Y ^Pn{S,i). (2) 

This index has certain appealing properties. In particular, Yli^x ^f 1 = PD(X), and 
there is a simple formula for quickly computing ifjf 1 . The index ifj sh also has a stochastic 
interpretation, but this is not based on extinction or survival of species, rather on the 
expected contribution to PD of each species under all possible orderings of the total set 
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of species (for details see Haake et al. 2005). The index ip sh allocates existing PD 'fairly' 
amongst the species, whereas tp quantifies the expected contribution of each species to 
future PD. 

3. Computing the HED index 

Computing the HED index directly via (pQ) could be problematic as it requires 
summation over all the subsets of X — {i} and this grows exponentially with |X|. However 
we now show that the index can be readily and quickly computed, both for rooted and 
unrooted trees. This polynomial-time algorithm for computing tjj thus complements (but 
is quite different to) the polynomial-time algorithm described by Haake et al. (2005) for 
computing -^ sh . 

3.1. Rooted trees 

For a rooted phylogenetic X-tree T let C(e) denote the set of species in X that are 
descended from e (i.e. the clade that results from deleting e from T). For i e X, let 
ei, 62, 6k (k = k(i) > 1) denote the edges (branches) on the path from i to the root of T, 
listed in the order they are visited by that path. Recall that /(e) denotes the length of edge 
e. The proof of the following theorem is given in the Appendix. 

Theorem 3.1 

k 

^ = X! i ^ n e J 

r=l jeC(e r )-{i} 

Note that in this (and the next) theorem we adopt the convention Ylje® e i = ^> w hich is 
relevant for the first term (r = 1) in the sum as C(e r ) — {i} is empty. Thus the first term 
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in the summation expression for ipi given by Theorem 13.11 is simply Z(ei), the length of the 
pendant edge of T incident with species i. 

3.2. Example 

We can apply the HED index to the members of the rooted tree depicted in Fig. 1. For 
example, to compute vpA by using Theorem 13.11 we have tpA — 1 + 1 • e.B + 1 • tB^c = 1.19. 
By inspection, we can see that the most valuable species will be D, since it shares an 
edge with only one other species above the root, and that this species (E) has a high 
P(ext). At the other extreme, A shares its path to the root with two other species, and 
one of them (C) has a low P(ext). It should therefore receive a low HED value. The 
computed values are ipu = 2.9, iJjb = 2.71, ipE = 2.1, ijjc = 2.09, and ipA = 1.19. Using 
the Shapley index (Haake et al. 2005), D and E are ranked first (with value = 2.63), 
followed by C (2.33) and then A and B (1.75). Pavoine's QE metric (Pavoine et al. 2005) 
returns the same ranking as does the Shapley. A portal for computing HED is available at 
http://www.disconti.nu/-phylo/emd.dpf 



3.3. Unrooted trees, and properties of the index 

We now provide a similar formula for efficiently computing the HED index for unrooted 
trees. Given a leaf i of T and an edge e of T, e induces a split of X into two disjoint 
subsets, and one of these subsets, which we denote as Cj(e), contains i. The proof of the 
following theorem is given in the Appendix. 
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Theorem 3.2 

e ^)-( n ^-c 1 - n ^ 

e&E{T) j e d(e)-{i} jeX-C^e) 

Notice that the rooted HED index is just a special case of the unrooted HED index 
(indeed Theorem 13. II can be deduced from Theorem 13.21) . To see this, given a rooted tree T 
attach a new leaf p to the root via a new edge to obtain an unrooted tree, and assign the 
new edge length 0. Let e p = 0. Then it is easily seen that the HED index for T is just the 
HED index for the derived unrooted tree. 

Using Theorem 13.11 it can be shown that if T is a rooted phylogenetic tree then the 
condition: 

J2^ = PD(X), (3) 

iex 

holds for all selections of positive branch lengths and e's if and only if T is a 'star tree' 
(that is, every leaf is adjacent to the root). Moreover Theorem 13.21 shows that there is 
no unrooted phylogenetic tree T for which (jHJ) holds for all positive branch lengths and e 
values (of course d3J) may hold on phylogenetic trees - either rooted or unrooted - if the 
branch lengths and e values take certain values). This contrasts with the index ip sh which 
satisfies J2iex ^f 1 = PD(X) on all unrooted phylogenetic trees and choices of branch 
lengths, a property that is referred to as the Pareto efficiency axiom by Haake et al. (2005). 
In the setting of this paper we should not be surprised that (j3J) holds for ip only in very 
special cases since we are not trying to divide out existing PD amongst present taxa (one 
motivation behind ip sh ) but rather quantify the expected contribution each species makes 
to future PD. 



4. Application 



We compared the HED index with the Shapley (Haake et al. 2005) values for the 
Prosimians (Mammalia: Primata), a group of approximately 50 species with a broad 
range of extinction probabilities. This group includes the Aye-Aye, the lemurs, the lorises 
and galagos. We made use of a recent dated Supertree of the order Primates (Vos and 
Mooers 2004; Vos 2006), see Fig. 2, and Red List risk designations from the IUCN 
(www.iucnredlist.org, accessed February 2006). Following Isaac et al. (2007) and Redding 
and Mooers (2006), we first converted the five categories of risk (CR, EN, VU, NT, and 
LC) to probabilities of extinction. Under the IUCN criteria, the species in the VU category 
are given a P(ext)=0.1 over the next 100 years. We gave the lowest and highest threat 
categories very conservative probabilities of extinction over the next 100 years of 0.001 and 
0.9 respectively, leaving P(ext) = 0.5 for EN, and P(ext) = 0.01 for NT: this scale is very 
similar to that calculated from real population viability analyses for birds (Redding and 
Mooers 2006). We are primarily interested in how the ranking of species changes using 
different approaches. 

The bivariate correlation between the metrics is high (0.94). Both measures chose 
the Aye-Aye (Daubentonia madagascarensis) as the most important species, followed 
by Perodicticus potto. Interestingly, the three most highly ranked species under current 
conservation policy (the critically endangered lemurs Propithecus tattersalli, Hapalemur 
simus, H. aureus) are nested well up in the tree (Figure 2) such that none of them were 
chosen in the top ten for either SV or HED. If we compare the rest of the rankings for these 
two metrics, the largest single difference is for the two Arctocebus species: they rank twelfth 
under SV (being relatively isolated on the tree), but only twenty-sixth under HED: because 
neither is severely threatened, the chances are good that their common path will persist. 
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Daubentonia madagascariensis EN 

— Phanev furcifer NT 
Cheirogaleus medius LC 

— Cheirogaleus major LC 

— Allocebus trichotis EN 

— Microcebus coquereli VU 
Microcebus murinus LC 
Microcebus rufus LC 

— lndri indri EN 
Avahi laniger NT 

— Propithecus diadema EN 
Propithecus tattersalli CR 
Propithecus verreauxi VU 

— Varecia variegata EN 
Eulemur coronatus VU 
Eulemur macaco VU 

— Eulemur fulvus LC 
Eulemur mongoz VU 

— Eulemur rubriventer VU 
Lemur catta VU 
Hapalemur simus CR 

— Hapalemur griseus LC 
Hapalemur aureus CR 
Lepilemur leucopus NT 

— Lepilemur mustelimis NT 
Lepilemur dorsalis VU 

— Lepilemur edwardsi NT 
Lepilemur microdon NT 
Lepilemur ruficaudatus NT 

— Lepilemur septentrionalis NT 
Perodicticus potto LC 

— Loris tardigradus EN 

— Nycticebus coucang LC 

— Nycticebus pygmaeus VU 
Arctocebus calabarensis NT 
Arctocebus aureus NT 
Galago senegalensis LC 

— Galago alleni NT 
Galagoides demidoff LC 

— Galagoides zanzibaricus NT 
Galago moholi LC 
Galago gallarum NT 
Galago matschiei NT 
Otolemur crassicaudatus LC 

— Otolemur garnettii LC 
Euoticus elegantulus NT 
Euoticus pallidus NT 

Fig. 2. — Prosimian species tree and associated IUCN threat categories. CR: critically 
endangered, P(ext) = 0.9; EN: endangered, P(ext) = 0.5; VU: vulnerable, P(ext) = 0.1; 
NT: near threatened, P(ext) = 0.01; LC: least concern, P(ext) = 0.001. Edge lengths are 
on an arbitrary scale that represents time since divergence. 
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Both measures are very heavily influenced by the pendant edge (PE) length of the 
focal species (with correlations PE vs. SV=0.94, and PE vs. HED=0.98). PE { is always 
part of the marginal increase to S, while interior nodes are most likely represented with 
high probability, especially for larger and more balanced trees. PE is, however, a poor 
predictor of HED for Propithecus, for H. griseus, and for Nycticebus couang (Figure 
2). The first two groups contain the three most endangered species, increasing the value 
of close relatives. Nycticebus is an isolated genus, and N.couang's sister species is listed 
as vulnerable ( P(ext) = 0.1). Likewise for Propithecusverreauxi - although it has close 
relatives and so a short PE, these relatives are at high risk of extinction, which increases its 
value; this is what we saw with species D in figure 1. 

4.1. Incorporating the focal taxon's extinction risk (HEDGE scores) 

The effect of close relatives' risk status on one's own value is precisely the strength of 
the HED approach. However, the fact that the extinction risk of other species affects a 
focal species, but its own risk does not is somewhat counter-intuitive. We address this by 
showing how it is possible to write the HED index as the sum of two terms each of which 
takes into account the extinction risk of the focal species. To describe this further, let I be 
the random variable which takes the value i if the focal species i survives (at the future 
time under consideration) and which otherwise takes the value of the emptyset (i.e. 0) if i 
goes extinct. 

Let 

$ = E[PD(Si U {i}) - PD(S t U /)], 

where, as before Si is the random subset of species in X — {i} that survive. In words, ^[ is 
the increase in the expected PD score if we condition on the event that species % survives. 
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Similarly, let 

^ = E[PD(S l Ul)-PD(S t )}. 

In words, ip" is the decrease in the expected PD score if we condition on the event that 
species % becomes extinct. The following result describes how to compute these two indices 
easily from the HED index, and verifies that they add together to give the HED index (its 
proof is given in the Appendix). 

Theorem 4.1 

(i) ip[ = e i -ijj i , 

(ii) ^ = (1 - Ci ) • Vi, 

(iii) IP'+IP'^IP,. 

The approach of assigning a value to a species which is a function of its phylogenetic 
distinctiveness and its extinction probability has been referred to as 'expected loss' by 
Redding and Mooers (2006) and, more evocatively, an 'EDGE' score (Evolutionarily 
Distinct and Globally Endangered) by Isaac et al. (2007). 

In the same spirit we will call ip' and ip" (which extend our HED index ip) HEDGE 
(heightened evolutionary distinctiveness and globally endangered) scores. The HEDGE score 
ip't is more relevant when evaluating actions that might save species, whereas the HEDGE 
score ip" is appropriate when evaluating actions that might cause the extinction of species 
(such as building a dam). 

One potential advantage of HED and HEDGE over previous scores is their flexibility 
in designing conservation scenarios. So for instance, we can choose IUCN-ranked species 
for which conservation is cheap and/or already partially successful, set their P(ext) to 0, 
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and see how rankings of other species change. Alternatively, we might want to increase the 
P(ext) to 1.0 for certain species to see how others are affected. 

Most generally, HED and HEDGE could be incorporated in an assessment of species 
value that included many factors besides risk and future contribution, e.g. the ecological, 
distributional and aesthetic values enumerated by Avise (2005), and the costs of recovery 
and probability of its success. 
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5. Appendix: Proofs of theorems. 

Proof of Theorem \3.1\ 

First observe that the only edge lengths that contribute to A PD (S,i) are those from the set 
{/(ei),/(e 2 ), . . . ,l(e k )}. 

Consequently, for the random set Si of surviving species of X — {i} we have 

A PD (Si,i) = J~] l{e r ) ■ I r (Si) 

r>l 

where I r {<Si) is the 0, 1 indicator random variable that takes the value 1 precisely if e r is 
not an edge of the subtree of T connecting the taxa in Si and the root of T; since this is 
the only situation for which e r lies in the subtree of T connecting Si U {i} but not in the 
subtree of T connecting Si. 

Thus, by linearity of expectation, 

fa = E[A PD (S h i)} = E[^/(e r ) • I r (Si)] = ^ l(e r ) ■ E[/ r («Sj)], 

r>l r>l 

and since I r (Si) takes the values and 1, E[/ r («Sj)] = F[I r (Si) = 1]. Thus, 

^ = ^2l{e r )-¥[I r {S i ) = l]. 

r>l 

Now, the event l I r (Si) = V occurs precisely if none of the elements in C(e r ) — {i} survive, 
and this latter event has probability njec(e r )-{i} e r Substituting this into the previous 
equation establishes the theorem. □ 

Proof of Theorem \3.S\ 

For % G X and the random subset Si C X — {i}, we have 

A PD {Si,i)= l (e)-h{Si) 
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where I e {<Si) is the 0, 1 indicator random variable taking the value 1 precisely if Si consists 
of no elements of Cj(e) — {i} and at least one element of X — Cj(e). 

Thus, 

e£E(T) 

and by the independence assumption 

F[I e (S t ) = 1] = P[5, n (die) - {t}) = 0] • F[Si fl (X - die)) + 0] 

and so 

F[usi) = i] = ( n ^-m 1 - n ^ 

i6Ci(e)-{i} ieX-QCe) 

as claimed. □ 



Proof of Theorem 4-1 



By definition ^ = E[PP>(<S; U {z}) - PD(Si)], and so E[PD(5i U {i})] = ^ + E[PD («%)]. 

Now we can write the unconditional expectation E[PD(Si U /)] as the weighted sum of 
conditional expectations E[PD(<S, U J)|J = {z}]^ = {*}) + E[PD(<Si U /)|J = 0]P(J = 0) 
and so 

E[PD(S, U /)] = (1 - ei )E[PD(5, U {i})] + eiE[PD(Si)]. 

Parts (i) and (ii) now follow by applying these equations (and the linearity of expectation) 
to the definitions of ip^ and Part (iii) follows directly from parts (i) and (ii). □ 



6. Acknowledgments 

This research was supported by a grant from the Marsden Fund, New Zealand to MS 
and AOM. AOM received additional support from NSERC Canada and the Institute for 



-17- 



Advanced Study, Berlin. We thank Dave Redding, Klaas Hartmann, and Andy Purvis for 
discussion. 



-18- 



REFERENCES 

Altschul, S. F., and D. J. Lipman. 1990. Equal animals. Nature 348:493-494. 

Avise J.C 2005. Phylogenetic units and currencies above and below the species level. Pages 
77-100 in Purvis A, Gittleman JL, and Brooks T, eds. Phylogeny and conservation 
Cambridge University Press, Cambridge, UK. 

Cousins, S.H. 1991. Species diversity measurement: choosing the right index. Trends in 
Ecology and Evolution 6:190-192. 

Faith, D.P. 1992. Conservation evaluation and phylogenetic diversity. Biological Conservation 
61:1-10. 

Haake, C-J, Kashiwada, A., and F.E. Su. 2005. The shapley value of phylogenetic trees. 



IMW Working Paper 363 (363). ArXiv (q-bio.QM/0506034) 



Hartmann, K. and M. Steel. 2007. Phylogenetic diversity: from combinatorics to ecology. 
Forthcoming. In O. Gascuel and M. Steel eds. Reconstruction evolution: New 
mathematical and computational advances. Oxford University Press. 

Isaac, N.J.B., Turvey, S.T., Collen, B., Waterman, C, and J.E.M. Baillie. 2007. Mammals 
on the EDGE: conservation priorities based on threat and phylogeny. PLoS ONE 
2(3): e296. 

May, R.M. 1990. Taxonomy as Destiny. Nature 347:129-130. 

Mooers, A.O. and S.B. Heard. 1997. Inferring process from phylogenetic tree shape. 
Quarterly Review of Biology 72:31-54. 

Myers, N.,Mittermeier, R. A., Mittermeier, C. G., da Fonseca, G. A. B., and J. Kent. 2000. 
Biodiversity hotspots for conservation priorities. Nature 403:853-858. 



-19- 



Pavoine, S., Oilier, S., and D. Pontier. 2005. Measuring diversity from dissimilarities with 
Rao's quadratic entropy: Are any dissimilarities suitable? Theoretical Population 
Biology 67:231-239. 

Pavoine, S., S. Oilier, and A.B. Dufour. 2005. Is the originality of a species measurable? 
Ecology Letters 8:579-586. 

Redding, D.W. and A.O. Mooers. 2006. Incorporating evolutionary measures into 
conservation prioritization. Conservation Biology 20:1670-1678. 

SARA. 2002. Bill C-5, An Act respecting the protection of wildlife species at risk in Canada. 
Government of Canada, Ottawa. 

Shapley, L.S. 1953. A value for n-person games. Annals of Mathematical Studies 28:307-317. 

Vos, R.A., and A.O. Mooers. 2004. Reconstructing divergence times for supertrees: a 

molecular approach. Pages 281-299 in O. R. P. Bininda-Emonds, ed. Phylogenetic 
supertrees: combining information to reveal the Tree of Life. Kluwer Academic Press, 
765. Dordrecht. 

Vos, R. 2006. Inferring large phylogenies: the big tree problem. PhD. Thesis, Simon Fraser 
University. 



This manuscript was prepared with the AAS IATjrpC macros v5.2. 



